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Abstract 

We apply the strategy of regions within dimensional regularization to find functions 
involved in evolution equations which govern the asymptotic dynamics of the Abelian 
form factor and four-fermion amplitude in the SU(N) gauge theory in the Sudakov 
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production in the e + e~ annihilation at high energy. 
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1 Introduction 



The asymptotic behavior of various amplitudes in the Sudakov limit has been investi- 
gated in QED and QCD, with summation of the leading double 0, |3|, [| and subleading 
H |, [7|, ||, [|, [10], |Tl| logarithms. Evolution equations that govern the dynamics of the am- 
plitudes in the Sudakov limit have been obtained in refs. || ^ || |10|| . In the present paper, 



we apply this approach to the next-to-leading analysis of the Abelian form factor and the 
four fermion amplitude in the SU(N) gauge theory. We evaluate functions that enter the 
evolution equations in the next-to-leading logarithmic approximation by using, as an input, 
asymptotic expansions of one-loop diagrams. Solving these equations we sum up the lead- 
ing and subleading Sudakov logarithms. The expansions of one-loop diagrams are obtained 
by use of the so-called generalized strategy of regions [0 (see also |13| ) which enables us 
to systematically identify the nature of various contributions and the origin of logarithms. 
This strategy is based on expanding integrands of Feynman integrals in typical regions and 
extending the integration domains to the whole space of the loop momenta so that a crucial 
difference with respect to the standard approach || [7|, |[ ||, [UJ is the absence of cut-offs that 
specify the regions in individual terms of the expansions. This approach is applied within 
dimensional regularization [|TJJ] when all integrals without scale (not only massless vacuum 
integrals) are put to zero. 

We apply our results for the subleading Sudakov logarithms to the analysis of the dom- 
inant electroweak corrections to the process of the fermion-antifermion pair production in 
the e + e~ annihilation. In the standard model of weak interactions the W and Z bosons get 
their masses via the Higgs mechanism and the Sudakov logarithms naturally appear in the 
virtual electroweak corrections |0|. They grow rapidly with energy and become dominant 
in the TeV region available at the LHC or the Next Linear Collider. The analysis of the 
Sudakov corrections is thus of high importance for the next generation of accelerators. The 
leading and subleading electroweak Sudakov logarithms were discussed in 1T6L K7\ 



m one 



loop approximation. The effect of higher order leading logarithms was estimated in 
by computing, in a physical gauge, a contribution related to the multiple virtual Z and W 
bosons exchanges. A complete analysis of the leading logarithms in exclusive and inclusive 
electroweak processes was done in [|Tj|] on the basis of the infrared evolution equation ap- 
proach. The subleading Sudakov logarithms, however, turn out to be significant and should 
be taken into account to get a reliable estimate of high order corrections. 

The paper is organized as follows. In the next section, the Abelian form factor is analyzed 
in two typical Sudakov type regimes. The analysis is then extended in Section 3 to the four- 
fermion amplitude. The electroweak Sudakov corrections are considered in Section 4. We 
present our conclusions in the last section. 
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2 The Abelian form factor in the Sudakov limit 



The (vector) form factor which determines the amplitude of the fermion scattering in the 
external Abelian field in the Born approximation can be written as follows 



B 



X 



where the four-vector index on the left hand side is suppressed, p± is incoming and p 2 is 
outgoing momentum. 



There are two "standard" regimes of the Sudakov limit s = (pi — p 2 ) 



-oo 



El: 



(i) On-shell massless fermions, p\ = p\ = 0, and gauge bosons with a small non-zero mass 
M 2 <C — s. Let us choose, for convenience, pi >2 = (Q/2, 0, 0, =F<5/2) so that 2pip 2 = Q 2 = — s. 

(ii) The massless gauge bosons and off-shell massless fermions p\ = p\ = — M 2 , M 2 <C — s. 
We choose p% = p\, 2 — (M 2 /Q 2 )p 2 ,i, where are defined as p\^ in the regime (i). 

The asymptotics of the form factor in the Sudakov limit can be found by solving the 
corresponding evolution equation 0. For the non- Abelian gauge theory, this equation was 
first derived in by factorizing collinear logarithms in the axial gauge. In the first regime, 
it reads 

rQ 2 Hr 

- 7 (a(a;)) + C(a(Q 2 ))+e(«(M 2 )) 







d\nQ 2 
Its solution is 

F = F (a(M 2 ))exp 



M 2 X 

Q 2 dx 

M 2 X 



^ 7 (a(x')) + C(«(x))+e(«(M 2 )) 

M 2 T 



(2) 



(3) 



A generalization of eq. (0) to the regime (ii) was found in |10 



d 



«91nQ 2 ~ 
Its solution is 



Q 



dx 



M 4 /Q 2 X 



7 (a(x)) + C(«(g 2 )) + C'(«(M 4 /Q 2 )) + e(«(M 2 )) 



(4) 



{rQ 2 dx \ r x dx 1 
— 9 —l(a(x')) + ((a(x)) + aa(M 2 )) 
JM 2 X [JM 2 X' 



+ 



M 2 



dx 



M 4 /Q 2 X 



M 2 dx > 



X 



(5) 



The functions Fq and ^ are, generally, different in the two regimes. We are interested in the 
next-to-leading logarithms. Therefore we should keep renormalization group corrections to 
the leading logarithmic approximation as well as single infrared and renormalization group 
logarithms. In this approximation, the form factor (i) takes the form 



F = F (a) exp 



Q 2 dx r x dx' 



M 2 X JM 2 X 



^ 7 («(x')) + (C(«) + ^(«))ln(Q 2 /M 2 ) 



(6) 
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and, in the regime (ii), we have 



„ / x I fQ 2 dx r x dx' , . ... r M ' 2 dx r M2 dx' , . ... 
F = F (a)exp / — / — 7 a x' + / — / — 7 a x' 

\Jm 2 X Jm 2 X' J M i /Q 2 X Jx X' 



+ (C(a) + C , («)+^(«))ln(Q 2 /M 2 ) . (7) 



All the functions in the exponent have to be computed in one loop, and the one loop running 
of a in the argument of the 7 function should be taken into account. Note that in the next- 
to-leading order we cannot separate the functions ( and £ but we will see that vanishes in 
one loop. 

In the covariant gauge, the self energy insertions to the external fermion lines do not give 
Q-dependent contributions. The one loop calculation of the vertex correction gives 

01 

F = tt C f (~ V o + 2V i + 2(1 - 2e)^ 2 - KQ F B , (8) 

Z7T 

where Cf = {N 2 — 1)/(2N) is the quadratic Casimir operator of the fundamental represen- 
tation of SU(N) group and the functions involved in are given by 

f7r d/2 e - 7E 6 s -iy o ^ 
l ^ e -^s~\p 1 +p 2 ) t y l , 

-I- ( n. <— » iA 

■ (9) 



d d k 


(k 2 


- 2 Pl k){k 2 - 2p 2 k){k 2 


— M 2 ) 




d d k k^ 




(k 2 


- 2 Pl k)(k 2 - 2p 2 k)(k 2 


— M 2 ) 








(k 2 


-2 Pl k){k 2 -2p 2 k)(k 2 


— M 2 ) 



i7l d/2 e - lE e 



9/ivV 2 H V 2 



in the regime (i). The corresponding relations for the regime (ii) are obtained from eq. (§) 
by the following substitution 

ik 2 - 2 Pl k){k 2 - 2p 2 k){k 2 - M 2 ) -> (k 2 - 2 Pl k - M 2 ){k 2 - 2p 2 k - M 2 )k 2 , (10) 

with another set of the functions V involved. We omit iO in k 2 — 2p\k + iO, etc. for brevity. 
Here 7e is the Euler constant. We work in dimensional regularization [14] with d = 4 — 2e. 
We also usually omit the factor (/i 2 ) e per loop and write it down only in the argument of the 
renormalization group logarithm. Note that only the tensor structures giving unsuppressed 
contributions to the form factor are kept in the representation of the third integral. 

To expand these integrals in the limit Q 2 3> M 2 we apply a generalized strategy of regions 
formulated in |T2[ and discussed using characteristic two- loop examples in [TBI]: 



Consider various regions of the loop momenta and expand, in every region, the inte- 
grand in Taylor series with respect to the parameters that are there considered small; 

Integrate the expanded integrand over the whole integration domain of the loop mo- 
menta; 
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Put to zero any scaleless integral. 



The following regions happen to be typical in the Sudakov limit [20 



hard (h) 

1- collinear (lc) 

2- collinear (2c) 

soft (s) 
ultrasoft (us) 



k~Q, 

k + ~ Q, A;_ ~ M 2 /Q , k ~ M , 
A;_ ~ Q, fc + ~ M 2 /Q , ~ M, 
fc ~ M, 
fc ~ M 2 /Q. 



Here fc± = /c ± &3, A; = (fei, A^)- We mean by ~ Q, etc. that any component of k^ is of 
order Q. 

Keeping the leading power in the expansion in the limit Q 2 /M 2 — > 00 we have in the 
regime (z)Q 



yh 



v c 



to 



1 1 ^ 
— In 1 



-In 2 



7T 

12 



1 1 1 civr 2 

+ _ l n (_ a ) _ l n (M 2 ) In (-s) + - In 2 (M 2 ) + — 
e z e 2 12 



2 ln [W 2 ) + 3 



Vi 



— +ln(-s) -2 

e 

--In (M 2 ) + l, 



(12) 



F 2 = tf* 



1 

2 ' 



ln 



(13) 



We denote by the index c the sum of the lc and 2c contributions. The pole in eq. (|1^) is of 
ultraviolet nature. It is not canceled by the collinear or ultrasoft contributions. Note that 
the corresponding ultraviolet renormalization group logarithm ln (— contributes only 
to the function (. 

1 In fact we do not need, in the LO and NLO, finite parts present in these and similar results. They would 
be, however, needed for the NNLO calculations. 
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The hard part in the regime (ii) is the same. The new ingredients read 



V C = -4 + -ln(M 2 )-ln 2 (M 2 ) + ^, 
V™ = ^ + -(ln(- S )-21n(M 2 )) + iln 2 (-s) 



e 2 £ ^ ' v 'J 2 

21n(M 2 )ln(-s) + 21n(M 2 ) + '— , (14) 



TV 2 
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Vn = In 2 



TV 2 



(m 2 ) ' 3 ' 



V? = - + 2 - In (M 2 ) , 
e 

v ' = ln S) ' (15) 

Soft regions generate only zero contributions (at least in the leading power). 
From the one-loop result we find 

7(«) = ~Cf£ ■ (16) 

Moreover it is clear from the above expressions that in the regime (i) the total double 
logarithms of Q come from the hard region while in the regime (ii) one half of the double 
logarithmic contribution comes from the ultrasoft region. This explicitly determines the scale 
of the coupling constant in the second order logarithmic derivative of the form factor in Q. 
It is Q in the regime (i) and M in the regime (ii). Furthermore, all the Q dependent terms 
in the ultrasoft contribution of eq. (H) are related to the 7 term and therefore ('(a) = 0. 



At the same time we cannot distinguish, in the one loop approximation, the contribution to 
the functions ( and £ coming from the collinear region because this region includes both Q 
and M scales. For the sum of these functions we find 

C(«)+£(«)=3C F £-. (17) 

To complete the Q-independent part of the one loop corrections to the form factor one has to 
include the fermion wave function renormalization determined by the self energy insertions 
to the external lines. In the regime (i) this brings the factor 




and in the regime (ii) this gives 




(19) 



6 



The ultraviolet poles of eqs. ( |T8| , [T^) cancel the ultraviolet pole in the first line of eq. (|T3| ) due 
to Ward identity and nonrenormalization property of the conserved vector current. Finally, 
in the NLO logarithmic approximation, we find the form factor (i) to be 



B 



+3cdn (Q 2 /M 2 



and the form factor (ii) to be 



exp 



2tt 



Q 2 dx r x dx' 



M 2 X JM 2 X 



—a(x') 



(20) 



B 



1 , 

F 2tt \ 3 
l '- dx /•'" ' dx' 

M 4 /Q 2 X 




exp 



2tt 



Q 2 dx f x dx' . .. 
— / — a{x ) 

M 2 X JM 2 X' 



rM 2 rM Ht' 11 

- / — / —a(x') + 3a In (Q 2 /M 2 ) 1 . 
Jm*/q 2 x L x< v ; ^ ' '\\ 



(21) 



The single logarithmic term in the exponents of eqs. (|20| , |2T|) has the following decomposition 



3 In (Q 2 /M 2 ) = 4 In (Q 2 /M 2 ) (IR) - In (Q 2 /M 



!{RG) 



(22) 



where we explicitly separate the infrared and renormalization groupQ logarithms which are 
related to V\ and V 2 integrals correspondingly. These logarithms are of essentially different 
nature (see a discussion below). Eqs. (f2(| pl|) are in agreement with the result of refs. |9|, |10 



1|) are in agreement with the result of refs. 
It is useful to distinguish the soft and collinear poles in e (resulting in logarithms) in the 
hard contribution. The collinear logarithms in a physical (Coulomb or axial) gauge originate 
only from the self energy insertions into the external fermion lines 0, |6], [7], |8|, 21], |22|, ^3 
and therefore are universal i.e. independent of specific processes. In contrast to the soft 
divergences (i.e. infrared divergences that are local in momentum space), the divergences of 
this type arise from the integration over angle variables. Consider, for example, the integral 
V\. The power counting tells us that there is no soft divergence in this integral (we have A; 5 
in the numerator and only fc 4 in the denominator). However, for non-zero k with k 2 ~ 
which is collinear to pi or p 2 (when pik ~ or p 2 k ~ 0), the integrand blows up. This 
follows from the fact that the factor l/{(k 2 — 2pk)k 2 ), with any p 2 = 0, generates collinear 
divergences. When integrating this factor in k , take residues in the upper half plane. For 
example, taking a residue at ko = —\k\ leads to an integral with l/(pk) = l/(p°\k\(l — cos9)) 
where 9 is the angle between the spatial components. Thus, for small 8, we have a divergent 
integration over angles dcos8/(l — cos#) ~ dO/6. The second residue generates a similar 
divergent behaviour — this can be seen by the change k — ► p — k. 

Within the method of regions the total divergence of the collinear region in general 
cancels both the soft and collinear poles of the hard part. Hence it is not straightforward 



2 Although this logarithm originates from the integration over the virtual momentum region between M 
and Q scales and does not depend on fj, we call it "renormalization group" one because it is directly related 
and can be read off the renormalization group properties of the Abelian vertex and the fermion wave function. 
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to separate the collinear logarithms. Let us note that the collinear divergences (but not the 
soft ones) can be regularized by introducing a finite fermion mass. The reason is that the 
light-like vector k cannot be collinear to the space-like vector pi in this case. This enables us 
to distinguish the soft and collinear poles. Introducing a small fermion mass m but keeping 
zero gauge boson mass we get the collinear poles in the hard part canceled by the poles of 
the contribution from the collinear region leaving the logarithms while the soft poles in the 
hard part are not canceled. Thus one can determine the origin of the poles in the hard part 
and therefore the origin of the logarithms. For example, in the hard part of the integral V\ 
the pole is of the collinear origin since it is canceled in the regime m 7^ 0, M = by 

= - - In (m 2 ) + 2 . (23) 

Thus the single infrared logarithm in eq. (|22]) is of collinear origin and therefore is universal. 
We should emphasize that this is not true for the renormalization group logarithm of this 
equation which depends on a specific amplitude and a model. For example it is different for 
the scalar form factor or for the vector form factor in a model with an additional Yukawa 
interaction of the fermions with the scalar bosons. 

A less trivial example is the integral Vq. In the hard part of Vq in this regime, the collinear 
part of the double pole is canceled by 

V C = ~ + ~ In (m 2 ) - ~ In 2 (m 2 ) - ^ (24) 
and transforms into the logarithm of s/m 2 but the soft single pole is left 

Vq = ~- In ft) + In ft) In (s) - \ In 2 ft) - £ . (25) 

3 The four fermion amplitude 

We study the limit of the fixed-angle scattering when all the invariant energy and momentum 
transfers of the process are much larger than the typical mass scale of internal particles 
\s\ ~ \t\ ~ \u\ ^> M 2 . Besides the extra kinematical variable the analysis of the four fermion 
amplitude is more complicated by the presence of different "color" and Lorentz structures. 
The Born amplitude, for example, can be expanded in the basis of color/chiral amplitudes 

Ab = % -^A" = % ^—Tp (-- (A d LL + A d LR ) + A C LL + A C LR + (L <-> i?)) , (26) 

where 

All = ^lI^il^lI^A, (27) 

A C LR = fah^lL^Rl^SR 



and so on. Here t a is the SU(N) generator, p\, p 3 are incoming and P2, Pa outgoing momenta 
so that t = (pi — P4) 2 and u = (p± + p 3 ) 2 — — (s + 1). For the moment we consider a parity 
conserving theory. Hence only two chiral amplitudes are independent, for example, LL and 
LF^\ Similarly only two color amplitudes are independent, for example, A and d. 

Let us first compute the one loop corrections in the regime (i). The total contribution of 
the vertex type diagrams is 

^ (C F (-V + 2V0 + ^V + . . .) A B , (28) 

where the ellipsis stands for the contribution without infrared logarithms. In the vertex 
correction involving the gauge boson selfcoupling, a contribution of the form ( p3f ) appears 
with m replaced by M and we have used the fact that the pole (logarithmic) term of Vf is 
the same both in eqs. (|23|) and (|12"D. 
The direct (uncrossed) box gives 



ig 2 a 

S 4-7T 



(5S + 3B 1 + AB 2 - 3B 3 - 6B 4 + 2B 5 ) 
Tr\ .. .. I) 



X { ~ if) ^^^^^^^ + ^27/^1^47^3 

- (3fi + B 1 -B 3 - 2B 4 + 2B 5 ) (29) 
x I [c F - - jpj 4>2t a 'J f ,'y5?pi4>4t a <y fl <y 5 tp 3 + ^-^27/^75^1 ^47^75^3 

where Tp — 1/2 is the index of the fundamental representation, 

B = -is 7(1,1,1,1,1) , B 1 = st J(l,2,l,l,l) , B 2 = st J{1,1, 2,1,1), 

B 3 = st J(l,3,l,l,2) , 54 = ^7(2,1,2,1,2), B 5 = s 2 J(l, 1, 2, 2, 2) , (30) 

and the functions J(ai, a-i, a 3 , 0,4, n) are proportional to scalar box integrals with shifted 
dimensions: 

J( ai , 02, a 3 , 04, U) = t ^--2-a :i -a 4 -l+d/2+n JJ^ _ ^| 

i 

r d d+n k 

X J (k 2 ) a i(k 2 - 2p x k - m 2 )° 2 (fc 2 - 2p 2 k - m 2 ) a *(k 2 - 2(p 1 - p A )k + t) a * ' ^ ' 

We apply the generalized strategy of regions to obtain a table of asymptotic expansions of 
the basic integrals in the leading power. Note that, besides the hard contribution, there 
are two groups of the lc, 2c and us-contributions corresponding to two choices of the loop 
momentum when it is considered to be the momentum flowing through one of the two gauge 



3 A translation to the SM with the different structures of the amplitudes for different chiralities will be 
considered in Section 4. 
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boson lines of the given box. Keeping terms with the leading and subleading logarithms we 
have 

B = B 5 = Bt - B 3 - 2£? 4 = , (32) 

fl*(a,f) = -l + Iln(-t) + iln 2 (- S )-ln(- S )ln(-t) , 
e z e 2 

B' 2 ( 3 ,t) = i-iln(-t)-ln(M 2 )ln(-t)-iln 2 (M 2 ) , (33) 

*<•.«) - Mm9+'»SM; 

and the direct box contribution reads 



sw B 2 <s,t)^C F -f)A^C F fA'). (34) 
The crossed box contribution can be obtained in the same way: 

£ZB, { .,«)((c,-%-%y + C*A<), (35) 

where Ca = N is the quadratic Casimir operator of the adjoint representation. The rest of the 
one loop logarithmic contributions from the vertex corrections and the self-energy insertions 
are of the renormalization group nature. In addition to the vertex and external fermion self- 
energy contributions considered in the previous section the renormalization group logarithms 
set the scale of g in the Born amplitude to be Q. 

The total one loop correction in the logarithmic approximation reads 



s 2tt 



(36) 



Note that the next-to-leading logarithms do not depend on chirality and are the same both 
for the LL and LR amplitudes. 

Now the collinear logarithms can be separated from the total one-loop correction. For 
each fermion- ant if ermion pair, they form the exponential factor found in the previous section 
(eq. ( POD ). This factor in addition incorporates the renormalization group logarithms which 
are not absorbed by changing the normalization scale of the gauge coupling. The rest of the 
single logarithms in eq. fl3Tf ) is of the soft nature. Let us denote by A the amplitude with 
the collinear logarithms factored out. It can be represented as a vector in the basis A x , A d 
and satisfies the following evolution equation f 



q^qsA = X(*(Q 2 ))A , (37) 
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where \ is the matrix of the "soft" anomalous dimensions. From eq. (|37| ) we find the elements 
of this matrix to be, in units of a/(47r), 



X\d 



-2C A \n [-)+a( Cf -^) Inf? 



N 



The solution of eq. 

A = A\(a(M 2 )) exp 



XdX — 

Xdd = . 
reads 
Q 2 dx 



4^1n^ 

JV \t 



"" ' 7 



(38) 



, — Xi(a(a;)) 



A°(a(M 2 ))exp 



Q 2 dx 

, — X2(a(a)) 

M 2 X 



(39) 



where Xi are eigenvalues of the x matrix and A® are Q-independent vectors. Note that in 
higher orders the matrices x f° r different values of Q do not commute and the solution is 
given by the path-ordered exponent ||. 

Equations fl53|) imply that only the hard parts contribute to eq. (|37|). This fixes the 



scale of a in this equation to be Q. By this reason the matrix ( 3q) is the same in the 
regime (ii). Hence in the next-to-leading logarithmic approximation the difference between 
the corrections to the four-quark amplitude in the regimes (i) and (ii) is the one between 
the factors © and flU). 

In the Abelian case, there are no different color amplitudes and there is only one anoma- 
lous dimension 



X 



4 In 



(40) 



4 Sudakov logarithms in electroweak processes 

We are interested in the process /'/' - * f f ■ In the Born approximation, its amplitude is of 
the following form 

<4* = v E (^ + 4^)47, (4i) 

S I, J=L, R K 4 J 

where 

41 = TilMjl.fj , (42) 

tw = tan6*vi/ with 9 W being the Weinberg angle and Tf (Yf) is the isospin (hypercharge) of 
the fermion which depends on the fermion chirality. 

To analyze the electroweak correction to the above process we use the approximation 
with the W and Z bosons of the same mass M and massless quarks and leptons. A fictitious 
photon mass A has to be introduced to regularize the infrared divergences. Let us consider 
first the equal mass case A = M, where we can work in terms of the fields of unbroken phase 
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and the result of Sects. 2 and 3 for the regime (?) can be directly applied by projecting on a 
relevant initial/final state. For each fermion-antifermion pair the factor ( ^0[) takes the form 



exp 

where 



T f (T f + l)+t 2 w ^pj (L(s)-3l(s)) 



(43) 



16tt 2 \M 2 

< 44 > 



and we neglect the running of the coupling constant in the integral in eq. fl20|) but fix the 
scale of the coupling g and tw9 in the double logarithmic contribution to be Q. The soft 
anomalous dimension for / and/or J = R is Abelian and, in units of g 2 / (16ti 2 ) , reads 

X = t 2 w Y f ,Y f \n(^j . (45) 

The matrix of the soft anomalous dimensions for / = J = L is a sum of the Abelian and 
non- Abelian parts 

Xxx = -4 In (^ + (t 2 w Y f ,Y f + 2) In 
3 , (u 



Xx d = ~ 4 ln u 

Xd x = 41n(y) , (46) 



Xdd = tw Y f' Y f ln 



The photon is however massless, and the corresponding infrared divergent contributions 
should be accompanied by the real soft photon radiation integrated to some resolution energy 
uj Tes to get an infrared safe cross section independent on an auxiliary photon mass. At the 
same time the massive gauge bosons are supposed to be detected as separate particles. In 
practice, the resolution energy is much less than the W (Z) boson mass so the soft photon 
emission is of the QED nature. This cancels the infrared singularities of the QED virtual 
correction. We should therefore separate the QED virtual correction from the complete result 
computed with the photon of some mass A and then evaluate the QED virtual corrections 
together with the real soft photon radiation effects with vanishing A. It is convenient to 
subtract the QED contribution computed with the photon of the mass M from the obtained 
result for the virtual corrections and then take the limit A — > for the sum of QED virtual 
and real photon contributions to the total amplitude. In the language of the approach of 
ref . [|TI| , this prescription means that we use the auxiliary photon mass A as a variable of the 



evolution equation below the scale M and the subtraction fixes a relevant initial condition 
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for this differential equation. This leads to a modification of the factor 
anomalous dimensions |47|) . 

The common factor for each fermion-antifermion pair becomes 



and the soft 



exp 



- (r f (T f + 1) + 4-2 _ s * wQ A ( L ( s ) - 3J(«)) 



where sw = sin 6w Then we have 

X = {t 2 w Y f ,Y f - As 2 w Q r Q f ) In (~ 



(47) 



(4£ 



and the matrix of the soft anomalous dimension for / = J = L is 

Xaa = 



-4 In (£) + (f w Y f ,Y } - ^ w Q f ,Q f + 2) In 



XXd 
XdX 
Xdd 



--^1 

I hl( - 



(49) 



(^y, - ^g^g,) in 



Now we can estimate the dominant one- and two-loop logarithmic corrections. The renor- 
malization group logarithms which are not included into eq. ( fTTD can be trivially taken into 
account by choosing the relevant scale of the coupling constants in the Born amplitude. At 
the same time, the remaining logarithmic corrections are of the main interest because they 
are supposed to dominate the (still unknown) total two-loop electroweak corrections. The 
one-loop leading and subleading logarithms can be directly obtained from eq. (|36|) . The 
corresponding corrections to the chiral amplitudes read 



T f (T f + l)+t 2 w ^l-s 2 w Q} + (f~n 



TpTj + 1 



YfiYf 



w ' 



-L(s)+3l(s)) 



+ 



-41n(-)+2m(-)+ln(-)4y / ,y / 



+ ln(-) t' w Y f Y f -^ w Q f ,Q f 



T f' T f + * 



w ' 



TfTf + -ln{ T )6 IL 8j L 



u 



(50) 



100 J M7 



where 5n = 1 for I = L and zero otherwise. The two-loop leading (infrared) logarithms 
are determined by the second order term of the expansion of the double (soft x collinear) 
logarithmic part of the collinear factors (47). The corresponding corrections to the chiral 
amplitudes are 



I (r f (T f + 1) + t^S- - s* w Q) + (/ ~ /')) L\s)Ag 



(51) 
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The two-loop next-to-leading logarithms are generated by the interference between the first 
order terms of the expansion of the double (soft x collinear) and single (soft+collinear+re- 
normalization group) logarithmic exponents. The corresponding corrections to the chiral 
amplitudes are of the following form 



T f (T f + l)+t 2 w ^l-s 2 w Q 2 f + (f 



f) 



T f (T f + l)+t 2 Jf-s 2 w Q 2 + (f 



f) 



+ 



u 



+ ln,| 



41nl-j+21n^ 

twY r Y,-ist v Q,,Q, 



+ In ( - ) i^YjYj 



T},T) + t 2 w 



TpTf 



t;,t; ■ i- Yf,Yf 

3 



In 



j w~ 
u 

7 



4 

SilSjl 



(52) 



f'Yf 



L(s)l(s)A f u f . 



With the expression for the chiral amplitudes at hand, we can compute the leading and sub- 
leading logarithmic corrections to the basic observables for e + e~ — > ff using standard for- 
mulae. Though the above approximation is not formally valid for small angles 9 < Mj yfs we 
can integrate the differential cross section over all angles to get a result with the logarithmic 
accuracy. Let us, for example, consider the total cross sections of the quark- antiquark/yU + /i~ 
production in the e + e~ annihilation. In the two loop approximation, the leading and next- 
to-leading Sudakov corrections to the cross sections read 



v/cr B (e + e — > QQ) 
a/a B (e + e~ -> qq) 
c/cr B (e + e~ — > ji + ji~ 



1 + 5.30 /(s) - 1.66 L(s) - 12Ml(s)L(s) + 1.92 L 2 (s) , 

1 + 20.54Z(s) - 2.17L(s) - 53.72 l(s)L(s) + 2.79 L 2 (s) , (53) 

1 + 10.09 Z(s) - 1.39 L(s) - 21.66/(s)L(s) + 1.41L 2 (s) , 



where Q — u,c,t, q — d, s, b, and we use Syy = 0.232 for numerical estimates. Numerically, 
L(s) = 0.07 (0.11) and l(s) = 0.014 (0.017) respectively for = 1 TeV and 2 TeV. Here 
M = M\y has been chosen for the infrared cutoff and g 2 /16ir 2 = 2.7 • 10~ 3 for the SU(2) 
coupling. For physical applications, the running of g and s\y, the W and Z boson mass 
difference and the top quark mass effects are important in the one loop approximation. For- 
tunately the first order corrections are known exactly beyond the logarithmic approximation 

n. 

Clearly, for energies at 1 and 2 TeV the two loop corrections are huge and amount up 
respectively to 5% and 7%. There is a cancellation between the leading and subleading 
logarithms and for the above energy interval the subleading contribution even exceeds the 
leading one. The higher order leading and next-to-leading corrections however do not exceed 
1% level. They can be in principle resummed using the formulae given above. 

For completeness we give a numerical estimate of corrections to the cross section asym- 
metries. In the case of the forward-backward asymmetry (the difference of the cross section 
averaged over forward and backward semispheres with respect to the electron beam direction 
divided by the total cross section) we get 



A FB /A B FB \ 



e e 



QQ) 



1 - 1.23 l(s) - 0.09 L{s) + 0.11 l{s)L(s) + 0.12 L 2 (s) 
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A FB /A^ B {e + e- -> qq) = 1 + 7.16 l(s) - 0.14L(s) - l.Ul(s)L(s) + 0.02 L 2 (s) , 
A FB /A| B (e + e" ->//" l >~) = 1 + 5.48ZO) -0.04L(s) - 6.39 /(s)L(s) + 0.27L 2 (s) . 

For the left-right asymmetry (the difference of the cross sections of the left and right particles 
production divided by the total cross section) we obtain 

A LR /Af R (e + e- -> QQ) = 1+ 10.17/(s) - 2.77L(s) + 6.38 l(s)L(s) - 0.91 L 2 (s) , 
A LR / Al R (e + e~ -> gg) = 1+ 11.66/(s)- 1.08L(s) + 15.27/(s)L(s) - 0.77 L 2 (s) , 
A M /Af R (e + e- -> = 1 + 118.07 l{s) - 13.74 L(s) - 1.13 Z(s)L(s) - 0.78 L 2 (s) . 

Let us compare our results with results of previous analyses. Our result for the one loop 
double logarithmic contribution is in agreement with fl6| . However the result for the one 



loop single infrared logarithmic contribution differs from flT7 |. The reason is that, in [17 
only the diagrams with heavy virtual bosons have been taken into account. There is an 
infrared safe contribution of the diagram with the virtual massless photon where the heavy 
boson mass serves as an infrared regulator that should be taken into account to get a com- 
plete (exponential) result. In one- loop approximation, this contribution comes from the box 
diagrams with the photon and Z boson running inside the loop. One of the two collinear 
regions of these diagrams (see Section 3) gives an infrared safe contribution that should be 
taken into account. The contribution from this diagram combined with real radiation has 
been obtained in analytical form in ||26|| . Note that neglecting contributions of this type 
leads also to the breakdown of the exponentiation of the double logarithms starting from the 
two-loop approximation |18| . 



Our result for the two- loop double logarithmic contribution is in agreement with |TJ|. On 
the other hand, the coefficients in front of the two-loop leading logarithms in eq. ([53]) with 
a few percent accuracy coincide with the result of |H| where the photon contributions were 



not considered. This is related to the fact that the virtual photon contribution not included 
to the result of [0 is suppressed by a small factor s 2 ^. 

Note that our analysis implies the resolution energy for the real photon emission to 
be smaller than the heavy boson mass. The above result for the cross section should be 
multiplied by the standard factor which takes into account the soft photon emission and 
the pure QED virtual corrections. This factor depends logarithmically on s, u res and the 
initial/final fermion masses but not on M z ,w- If the resolution energy exceeds M z> w the 
analysis is more complicated due to the fact that the radiation of real photons is not of 
Poisson type because of its non-Abelian SU(2)l component. A complete analysis of this 
problem in the double logarithmic approximation is given in ref. |TjJ. However, as we have 
pointed out the effects of the non-Abelian component of the photon are numerically rather 
small. 

5 Conclusion 

In the present paper, we have analyzed the asymptotic behavior of the Abelian form factor 
and four fermion amplitude in the SU(N) gauge theory in two standard variants of the 
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Sudakov limit: with on-shell massless fermions and massive gauge bosons and with off- 
shell massless fermions and massless gauge bosons. The generalized strategy of regions 
and dimensional regularization were used to obtain the asymptotic expansions of one-loop 
diagrams that determine the structure of the evolution equations for the amplitudes in the 
Sudakov limit up to the next-to-leading logarithmic approximation. By integrating these 
equations the next-to-leading logarithms were summed up. The method can be directly 
extended to the next-to-next-to-leading logarithms. To do this one can apply expansions of 
two-loop Feynman integrals within the strategy of regions (in the case of the form factor see 



examples of expansions of master scalar integrals in [jO^, [27]]) and insert two- loop information 
into the evolution equations. 

We have applied our results to the analysis of the electroweak corrections to the pro- 
cess of fermion-antifermion pair production in the e + e~ annihilation. The two-loop leading 
and next-to-leading logarithmic corrections to the chiral amplitudes which are supposed to 
saturate the total two-loop electroweak corrections in TeV region have been obtained. The 
corresponding corrections to the total cross sections and asymmetries of the quark-antiquark 
and n~ production in the e + e~ annihilation have been found to be of a few percent mag- 
nitude at the energy of 1-2 TeV. The next-to-leading infrared logarithms are comparable and 
even exceed the leading ones at this scale. 
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